Detection of spin reversal and nutations through current measurements 
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The dynamics of a single spin embedded in a the tunnel junction between ferromagnetic contacts 
is strongly affected by the exchange coupling to the tunneling electrons. Moment reversal of the 
local spin induced by the bias voltage across the junction is shown to have a measurable effect on 
the tunneling current. Furthermore, the frequency of a harmonic bias voltage is picked up by the 
local spin dynamics and transferred back to the current generating a double frequency component. 
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Detection and manipulation of single spins is an im- 
portant field of science since it pushes the limits of quan- 
tum measurements. Single spins would also be objects 
for qubits and, thus, be crucial for quantum informa- 
tion technology. Potentially, spintronics will replace con- 
ventional electronics devices with spin analogues where 
manipulation, control, and read-out of spins enable func- 
tionality with no or little charge dynamics.^^ One main 
issue concerning measurements on the quantum limit is 
whether the local dynamics of quantum object generate 
measurable responses to externally applied fields. 

Current-induced magnetic s witching has been ad- 
dressed for macroscopic planaiP^BEEl and magneto- 
mechanicaP systems, and spin-dependent tunneling 
in systems with non-coUinear configured ferromag- 
netic leadSj'^'^Iini as well a s questions concerning 
spin-tran sfer torqu e j ^^*^^*^*^ * ^^ ! decoherence,'^ and spin 
nutationgl^IISIilEol microscopic spin systems. It was 
recently demonstrated that the nutations of a local spin 
embedded in the tunnel junction between ferromagnetic 
leads can be electrically controlled and undergo stimu- 
lated spin reversal, using short bias voltage pulses.*^ 

The purpose of this paper is to discuss the efi^ect of 
the local spin dynamics on the current flow through the 
system. In an arrangement with one ferromagnetic and 
one non-magnetic lead, the local spin tends to align with 
the magnetic moment of the ferromagnetic leads when 
this lead acts as the source, see Fig. [T] upper panel. This 
alignment leads to an enhanced tunneling through the 
junction since the local spin can mediate the tunneling 
electrons at a higher rate. If the ferromagnetic lead acts 
as the drain, however, the local spin align anti-parallel 
to the magnetic moment of the lead which decreases the 
tunneling through the junction. The local spin motion 
generates significant variations in the current. Second, 
we discuss the effect of a harmonic bias on systems ar- 
ranged with two ferromagnetic leads in a non-parallel 
configuration, see Fig. [l] lower panel. The frequency 
of the bias voltage is picked up by the local spin mo- 
tion and is transferred back to the current, generating 
current components with frequencies loq and 2u)q. The 
doubled frequency component is a fingerprint of the in- 
teraction between the tunneling electrons and the local 
spin in presence of the bias induced magnetic field. 

The back action of the tunnel current on the local spin 



and the resulting spin dynamics was studied in Ref. 1191 
Here, we implement the spin dynamics of the localized 
spin into the current and explore the prospect of detect- 
ing the spin dynamics in the current. Consider a single 
spin embedded in a tunnel junction between two ferro- 
magnetic leads, see Fig. [T] The Hamiltonian for the 
system can be written as 

H = Hl + Hr + Hs + Ht, (1) 

where = Ep(9)a ep(9)ac|,(,)^Cp(,)^, where 

'^p{q)a- creates (annihilates) an electron in 

the left (right) lead at the energy ep(g)CT and spin a =ti i- 
The spin Hamiltonian is simply Hs = — ffMsBo • S, where 
g and /is is the gyromagnetic ratio and Bohr magneton, 
respectively, whereas Bp is the external magnetic field 
acting on the local spin S. Finally, the tunneling 
Hamiltonian Ht = Y.pqap{4>aTapCqi3 + H.c), where 
a, 13 are spin indices whereas the tunneling operator 
Taf} = ToSa/s -I- TiS • (Taf3, with the tunneling rate Tq 
between the leads, the tunneling rate Ti between the 
leads and the spin, and the vector of Pauli spin matrices 

Whenever there is a voltage applied across the 




FIG. 1: (Color online) Schematic views of two possible set- 
ups relevant to the present study. In the upper panel, the 
local spin is located on a substrate surface for SP-STM mea- 
surement. The lower panel shows the local spin embedded in 
the tunnel junction between two ferromagnetic leads. 
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junction, the tunneling current is given by I(t) 
-ed/dt(NL) = ie{[NL,'HT]), which results in {h = 1) 



where 



I{t) = 2eRe^ / e 

a/3 



iS(t>(t,t') 



where the operator Aafiit) — J2pq Cpa(i)cg/3(t). Here, the 
transformed electronic annihilation operators 



(3) 



with Kl(^ji) = HL(i?) + (A*L(fl) + VFL(fl)(t))A^L(fl) and 
a^p{q)'^- unequal chemical po- 
tentials of the system leads to the bias voltage + 
WL{t) — ^'R—WR{t) = e[Vdc+Vac(i)], and we have defined 
m = e t jVA, + VUt')]dt' , and 5<P{t,t') = 
Here we will consider biasing of the system with 



6l(t,0=e j {Vdc + Vac[Q{T-To)-Q{T-T{)])dT, (4) 



with tq < Ti, and 



Vac COS WoT)dr, 



(5) 



For later reference, we also define ^p(q)a = ^p(q)cr ~A'l(_r)- 
Before deriving the expression for the current we 
briefly discuss the dynamics of the local spin, for which 
an equation was derived in Ref. 19 and found to be we ll 
described by the Landau-Lifshitz-Gilbert equatiorPS^U 

dS{t) _ ^(^^cfS(t) ^ g^^^ _^ g^sS{t) X B{t), (6) 



dt 



dt 



where B = Bq + b|j^^ + B|j^^ is the efi^ective magnetic 
field. The energies associated with the spin dynamics, 
Hlol 1 ^eV, are much smaller than the energies for the 
electronic processes which are of the order of at least 1 
meV. The damping factor a{t) ~ T^/D at low tempej]- 
atures, where 2D is the band widths in of the leads,™ 
and it is reasonable to believe that this is the case also 
for finite temperatures. The leads are assumed to have 
large D, hence, a becomes negligible and as a result Eq. 
^ can be analytically solved. 

The induced magnetic fields can be expressed aJ^^HH 



pqa 



xlm / e^[(«p<'-«.")(*-*')+5</'(t,t')](|i'z (7) 

' — OO 



g^iBB^^l{t) = S f {x(2,)(i,i')Kx-n,y) 

J — OQ 

+ [Ki'J{t,t') - Ki'J{t,t')]n,i]dt' (8) 



pqa 

y^gil{ip<.-igi>){t-t')+S<P{t,t')] 



(9) 



and 



pqa 

-[/(^Pff) - /(e,.)]e*(«---«-)(*-*')}e^'*(*'*') (10) 

The local spin-dynamics can be described by 
its corresponding classical motion, which is con- 
sistent with the parametrization S = Sn — 
S'(cosiy9sin0,sin(/jsin0,cos0). Assuming time- 

independent external magnetic field Bq — BqZ, 
results in the equations of motion 



dip 
~di 



d0 
dt 



+S I [Ki^J{t,t')-Kif{t,t')]dt' cos 



S / Ki^J{t,t')dt' sine. 



(11) 



The expression for the current is conveniently divided 
into three components, e.g. I{t) = Io{t) + Iz{t) + I±{t), 
and using that the local spin is semi-classically described 
these contributions can be written as 

Io{t) = 2eRe^ /* To'[/(Cp.) " /fe^)] 

J —OO 



pqa 



(12) 



which accounts for the direct tunneling between the 
leads, while 



h{t) = 4eRe^ / TqTi 

pqa 



2To 



cos e{t) 



cose{t) 



xifi^pa) - /(e,.)]e^(«-^«-)(*-*')+^'^(*'*')di', (13) 
is influenced by the local spin moment, and 

Ij^it) = 2eRe^ f sin^ 0(t) [/(^p,) - /(C,^)] 

pqa 

includes spin-flip of the tunneling electron caused by the 
presence of the local spin. The resulting current is clearly 
modulated by the presence of the spin through the varia- 
tion of the spin polar angle 0{t), and through its doubled 
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frequency 26{t). This property will be exploited later in 
this paper, and it will be shown that the nutations of the 
local spin which are induced by the bias voltage can be 
measured in the current as higher order harmonics. The 
expressions for the current also show that the azimuthal 
motion, (/?(<), has no direct influence on the current. This 
would also be expected for a spin influenced only by an 
external magnetic fleld Bq = B^i. Thus, we need only 
to determine the polar angle motion of the local spin, for 
a complete description of its effect on the current. The 
general solution for the polar angle motion can be written 

al^-\ o t- ft- ^0 S J* J*' Kil\t' ,t")dt"dt'\ 

= zarctani tan — e •'*o''-°° y ' ' 

where is the initial polar angle at the initial time to. 
For the bias 64>i{t,t'), Eq. Q, we have 



de{t) 
dt 



= -2T:STf cjlaNLaNR-AeVdMro - t) 
+ eVd, cos \eVac{t - ro)]e-(*-^«)/^ 



-e(i - n)] + 



[e(i ~ To) 



eV^ZcCOS [eKc(Ti - tq)] 
COS \eVarXt - Ti)] + eVdc ( 1 



Q{t-Ti)\smB(t), (16) 



-(t-ri)/i 



where the time-scale r relates to the electronic tunneling 
processes. Such processes are of the order 1 fs, which 
is much faster than the characteristic time-scale for the 
spin dynamics for biases ranging in 1 — 100 mV. For e.g. 
a 1 ns bias pulse of 1 mV, the critical time scale is ~ 100 
fs, which should be within the realms of present experi- 
mental state-of-the-art-technology for nanoscale systems. 
Physically, the time-scale r means that the induced mag- 
netic fleld is a retarded response to the bias across the 



junction. In Eq. (16 1, N, 



^L(R)cr 



is the spin a density of 



states (DOS) in theleft (right) lead. 

The characteristic time-scale for the polar angle mo- 
tion can be written as 1/tc ~ 27reFTf cr^^A^LCT-^flg-- 
Parametrizing the spin-polarized DOS -/VL(fl')(j = iVo(l + 
o'PL(i?))/2,^'' where -1 < Vl{K) < 1, assuming TiiVo 
0.1, PL — -~Pr — 1/2, and V ^ I mV, gives Tc ~ 5 ps, 
which is sufficiently short to switch the spin from being 
I to 1 within 1 ns. 

Assuming initially that Oq = t:, that is, the local spin 
moment is |, see Fig. |2]a), showing the time-dependence 
of the polar angle for the bias pulse corresponding to 
6(j)i{t,t'). Also, assume that V^c > 0, and that pl = 1/2, 
Pn = 0, for TiA^o ~ 0.1. The assumption is thus that one 




FIG. 2: (Color online) a) Time-dependence ol the polar an- 
gle 9{t) (solid) for the pulsed bias voltage V{t) = eVdc + 
eVaclO{t - ro) - e(t - n)], (ro, n) = (1, 2) ns (dotted) (zero 
bias indicated through dashed line), b) Time- dependence 
of the total current I{t) (bold), and the contributions Io{t) 
(faint), Iz{t) (dashed), and /±(t) (dotted). 



lead is ferromagnetic whereas the other is non-magnetic, 
which would be applicable to scanning tunneling mi- 
croscopy (STM) measurements with a spin-polarized tip. 
The contribution from 1^ is then flnite whereas I± ~ 0, 
see Fig. |2]b) showing the total current I{t) (bold), and 
the contributions /o(i) (faint), Iz{t) (dashed), and I±{t) 
(dotted). The influence of the stationary fleld tends to 
flip the spin moment to which can be understood 
from Eq. (15 1 where the exponent has a negative sign. 
The contribution from 



increases at the instant of the 
spin-flip, since the local spin then aligns with the major- 
ity spin in the source (left lead) and can thus mediate the 
tunneling electrons at a higher rate. Momentarily, as the 
local spin flips from J, to t, the polar angle 9{t) ^ 0, ±7r 
which implies that the contribution from Ij_ peaks during 
a short time period. The peak of this contribution will, 
however, be hidden by the sharp increase of Iz during 
the same time period. Nonetheless, this analysis shows 
that the total current through the system increases as a 
response to the spin-flip of the local spin, see Fig. [2]b) 
(bold). 

Now, assume that the amplitude of the bias pulse 
Vac < —Vdc < 0. This bias pulse causes the current to 
flow from the right to left, and the induced magnetic fleld, 
which affects the polar angle motion of the spin, works 
as to reverse the moment of local spin. This can be seen 
in the sign of the exponent in Eq. (15 1, which is posi- 



tive and, thus, signifles that the polar angle approaches 
TT. Thus, we assume that the bias pulse is sufficiently 
strong and/or long to stimulate the spin reversal (that is 
going from f to |), see Fig. [2] a). At the instant of the 
spin-flip, the contribution from Iz decreases since the lo- 
cal spin moment becomes opposite to spin-polarization 
of the drain (left lead), see Fig. 2]b) (dashed). The lo- 
cal spin mediates the tunneling electron at a slower rate. 
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FIG. 3: (Color online) a) Total current I{t) through the sys- 
tem biased with V{t) = Vdc + Vac cos wo^, ^^o^ = 5h ■ lO^/vr, 
and e) its Fourier transform showing a main peak at a; = ojq 
and a second peak at = 2uJo- b) The current contribution 
Io{t), and f) its Fourier transform showing a single peak at 
LJ = tJo- c), d) The current contributions Iz{t) and I±{t), re- 
spectively, and g), h) their corresponding Fourier transforms, 
both showing peaks a.t u> = luq, 2u)o- Panels e) - h) display 
the relative intensity of the contributions in arbitrary units. 
Here, the spin-polarization pl = ~Pr = 0.9, T\/To — 1/2, 
iYo^i ~ 0.1, eVdc ~ lOOOojo, and eVac ~ 5loo. 



Analogously as in the above discussion, the contribution 
from /_L peaks during the short time period of the spin 
reversal, see Fig. |2]b) (dotted), nevertheless, the net ef- 
fect of the spin-flip is a reduced amplitude of the current, 
see Fig. [2]b) (bold). As the pulse terminates, the bias 
again becomes positive which reverses the local spin mo- 
ment from J, to "f . It should also be noticed that while the 
contribution from /o(i) follows the biasing of the system, 
it is unaffected by the local spin dynamics, as expected, 
see Fig. [2]b) (faint). 

For a harmonic bias (502 (i,i'), the current response to 
the local spin displays higher order harmonics which are 
only present when the leads are ferromagnetic and in a 
non-parallel configuration. The harmonic modulation, 
with frequency 0{t), of the bias is picked up by the local 
spin, however, its presence is transferred to the current 
such that a second harmonic, 29 (t), appears in the cur- 
rent modulation. This can be directly seen in the ex- 
pressions for the currents Iz and I± , both of which have 
terms with the factors cos26'(i). 

The transverse induced magnetic field j'^^Ki^Jdt', 
and the currents Iz and I± all have the time- dependence 



^i{x+eVdc){t—t' )+ieVac{sin uigt—sin uiot' ) / ujo fj^^' r:^ 

oc 

Jn{eVac/(^Q)Jm{eVac/iOo)S{x + eVdc + nujo) 




FIG. 4: (Color online) a) Polar angle motion (solid) and 
cos6(t) (dotted) for the bias V{t) = Vdc + Vac cosujot, and b) 
spin motion. In panel b), an external magnetic field Bo = 1 T 
has been added in order to speed up the azimuthal rotation, 
whereas other parameters are as in Fig. [3] 



It should be noticed that although this time-dependence 
arise in the current contribution Io{t), the net effect of 
the bias still yields a single harmonic variation of this 
contribution, see Fig. [sjb) and f), illustrating the time- 
dependence of /q (t) and its Fourier transform. The plots 
clearly shows that the direct tunneling current acquires 
only the frequency of the bias and no higher harmonics. 
This can be understood from the fact that Io{t) is a direct 
kinetic response to the electric field and does not reflect 
any of the local (spin) dynamics in the tunnel junction. 

The polar angle motion resulting from the time- 
dependence in Eq. ( 17 1 is given by 



9{t) — 2 arctan ( tan -— JJ^ exp 



27TSTf 



y^^NLaNRs(Tl„Jn{eVarJujQ)Jra{e.Vac/^^o) 
a 

sinwo('^ ~ '™)^ 



LOQ(n — m) 



(18) 



xe 



zcjQ {rt—ra)t 



(17) 



showing that the nutations of the spin caused by the bias 
induced magnetic fleld are directly modulated by the har- 
monic bias. Moreover, from this expression of the polar 
angle motion, it is clear that the lower the frequency, 
i.e. Wo ^ eVdc/W, the larger the amplitude of the an- 
gular displacement. Therefore, low frequencies provide a 
greater response to the current since it is influenced by 
the spin nutations through sm9{t) and cos9{t). 

In order to acquire a large response from the spin nu- 
tations it is also required that the magnetic moments 
of the leads are non-parallel. Therefore, assume that 
Pi^ = ~Pii = 0.9. Moreover, in order to further am- 
plify the response from the spin nutations, we assume 
that Ti/To = 1/2, NqTi - 0.1, and let eVdc IOOOwq 
and eVac ^ 5ujq- Under those conditions, which are rea- 
sonable from an experimental view point, the amplitude 
of the polar angle {maxd ~ min0)/7r 0.05, Fig. |4] a) 
(solid), which is a sufficiently large displacement to gen- 
erate a significant contribution from cos6'(i), sin0(t), see 
Fig. |4]a) (dotted), which is transferred into the current 
contributions /z(t) and I±(t). 
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FIG. 5: Total a) time-dependent current I{t), and b) its 
Fourier transform !/(<jj)1 in the case of non-spin-polarized 
leads, e.g. Pl = Pr = 0. Other parameters are as in Fig. 

El 

The current contributions Iz{t) and I±{t), and their 
corresponding Fourier transforms, are plotted in Fig. [3] 
c), d), g), and h), which illustrate that their contribu- 
tions modify the shape of the total current c.f. panels 
a), c), and d) in Fig. [S] and that this modification arise 
from current contributions with frequency 2^0 5 c.f. pan- 
els e), g), and h). The contributions with this double 
frequency is a direct fingerprint of the (exchange) inter- 
action between tunneling electrons and the local spin, 
in combination with the spin-polarization of the leads. 
In case of non-spin-polarized leads there is no magnetic 
field induced from the biasing of the systerrP^l and, thus, 
the time-dependent modulation of the bias is not trans- 
ferred to the local spin dynamics. Under those circum- 
stances the current consist of a single frequency contri- 
bution, see Fig. [5] which illustrates the time-dependence 



a) and Fourier transform b) of the total current, which 
is understood since then the current is simply a kinetic 
response to the electric field with no influence from the 
local (spin) dynamics. 

In summary, the effect of a local spin, embedded in the 
tunnel junction between ferromagnetic leads, on the cur- 
rent through the system has been studied. It has been 
shown that the spin reversal of the local spin, stimulated 
by the bias voltage induced magnetic field in the junction, 
generates a measurable variation of the current. This fea- 
ture thus enable direct recording of the local spin reversal 
in the tunneling current. It has, moreover, been demon- 
strated that the spin nutations induced by a harmonic 
bias voltage with frequency loq is transferred back into 
the current, generating components with frequencies loq 
and 2ujq. The double frequency component is a direct 
response of the spin motion in the current, which van- 
ishes whenever the transversal induced magnetic field is 
absent. It was also shown that the 2loq component is suf- 
ficiently large to be measured within the realms of todays 
state-of-the-art nanotechnology. Experiments of the dis- 
cussed spin dynamics measurements would be extremely 
intriguing and useful for further advances within basic 
science of nanoscale physics. 
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the International Ten Bar Cafe. 
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